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ABSTRACT 

Millisecond x-ray pulsars have weak magnetic dipole moments of ~ 10 16 Tm 3 compared to 
ordinary X-ray pulsars with dipole moments of 10 2() T m 3 . For this reason a surrounding accre- 
tion disc can extend closer to the neutron star, and thus reach a higher temperature, at which 
the opacity is dominated by electron scattering and radiation pressure is strong. We com- 
pute the self-similar structure of such a geometrically thin axisymmetric accretion disc with 
an internal dynamo. Such models produce significantly stronger torques on the neutron star 
than models without dynamos, and can explain the strong spin variations in some millisecond 
X-ray pulsars. 

Key words: accretion, accretion discs - magnetic fields - MHD -stars: magnetic fields - 
stars: neutron - X-rays:binaries - pulsars: general 



1 INTRODUCTION 

The first accretion-powered milli second pulsar, SAX J18 08.4- 
3658, was discovered a decade ago iWiinands & van derKlisI 19981 : 
IChakrabartv & Morgan! 1 19981) . Since then several other millisec- 
ond X-ray pulsars with spin periods between 1.7 and 5.4 ms have 
been discovered (e.g. XTE J0 929-314, XTE J1 751-305, XTE 1807- 
294 and XTE J1814-338; see lWiinands1l2004l for a review). These 
neutron stars have considerably weaker magnetic dipole moments, 
~ 10 16 Tm 3 , than the conventional X-ray pulsars. The accretion 
disc then extends much closer to the neutron star, where it reaches 
such a high temperature that the opacity is dominated by elec- 
tron scattering and radiation pressure can grow stronger than the 
gas pressure. On the other hand the accretion rates in these sys- 
tems can be so low that the Alfven radius might be located out- 
side of the co-rotation radius of the neutron star. Conventionally 
it has been thought that such systems will eject the mass that is 
transf erred from the companion s tar through the propeller mecha- 
nism (111 arionov & Sunvaev 1975), th ough it has been poin t ed ou t 
bv lSpruit & Taarrj d 19931) and Rappapor t. Fregeau & Spruitl J2004I) 
that the velocity excess righ t outside the co-rotatio n radius is insuf- 
ficient to eject all the matter. Ispruit & Taaml ll993) found that when 
the Alfven radius is only slightly larger than the co-rotation radius 
then the disc will cycle between one state in which the matter is ac- 
cumulated outside the co-rotation radius an d another state in which 
it is a ccreted onto the star. Later iRappaport. Fregeau & Spruitl 
d2004h constructed a steady accretion disc model, in which the en- 
tire disc is located outside of the co-rotation radius. The magnetic 
field of the neutron star is penetrating this disc and is mediating an 
exchange of angular momentum between the acc r etion d isc and the 
neutron star as in the model bv lGhosh&Lambl jl979h . Since the 
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entire accretion disc is rotating more slowly than the neutron star 
in this case, the disc brakes the rotation of the neutron star. 

The models of Ghosh & Lambl ( 119791) and 
Rappa port. Fregeau & Spruitl d2004l) have ignored the possi- 
bility of an internal dynamo in the accretion disc, and have 
considered only the magnetic field that is generated in the disc 
by the mismatch in the rotation of the neutron star and the disc. 
However the magnetic field due to a dynamo can strengthen the 
coupli ng between the accretion disc and the neutron star, and 
we dTessema & Torkelsson! |2010l) have shown that the angular 
momentum exchange between the accretion disc and the neutron 
star can increase by more t han an order of magni tude compared 
to the standard model by I Ghosh & Lambl dl979l) . Furthermore 
the dynamo is insensitive to the rotation of the neutron star, and 
can therefore in principle provide a spin-up torque even at a low 
accretion rate. 

A complementary approach to study the magnetic interaction 
between an ac cretion disc and a s tar is represented by the numerical 
simulations bv lLong et alj j2008). In these simulations it is possible 
to study more complicated magnetic field geom etries than are ac - 
cessible to semi-analytic studies (see for instance lLong et ail 2008). 
but the resolution is so low that it is not possible to study the appear- 
ance of small-scale turbulence and the resulting turbulent transport 
and dynamo action in the simulation. Thus these simulations do not 
really address the viability of a dynamo in the disc. However there 
are other things one can learn from these simulations. They show 
how the matter is streaming via funnel flows from the accretion disc 
onto the magnetic poles of the neutron star. Furthermore they show 
that angular momentum is also lost from the accre tion disc through 
the jet s that are accelerated away from the disc. iRomanova et al.l 
(2009) make a distinction between a conical jet that transport most 
of the mass, and a less dense axial jet inside the conical jet. The 
latter component can dominate the energetics if the star is in the 
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propeller regime. It is also noteworthy that these simulations con- 
firm that some matter might still be accreted when the star is in the 
propeller regime. 

The purpose of this paper is to extend our previous study (Pa- 
per I) of an accretion disc with an internal dynamo around a magne- 
tised star to the discs around millisecond pulsars. The difference is 
that due to the weak magnetic field the disc extends to much smaller 
radii, at which it becomes so hot that a significant part of the opacity 
is due to electron scattering, and radiation pres sure might dominate 
over the gas pressure. Therefore in the spirit o flShakura & Sunvaevl 
dl973ll 1976b we divide the disc into an outer region with gas pres- 
sure and free-free opacity, a middle region with gas pressure and 
electron scattering opacity, and an inner region with radiation pres- 
sure and electron scattering. 

In Sect. 2 we derive a single ordinary differential equation, 
which is generally applicable for the radial structure of the accre- 
tion disc independently of the equation of state and the form of 
the opacity. We then introduce the equations of state and the opac- 
ity prescriptions for the different regions of the disc in Sect. 3 and 
present the numerical solutions in Sect. 4. The properties of these 
solutions and their application to the millisecond X-ray pulsars are 
discussed in Sect. 5, and we finally summarise our results in Sect. 
6. 



2 GOVERNING EQUATIONS 

2.1 The magnetohydrodynamic angular momentum balance 

Our approac h is an extension of the st andard model for a thin ac- 
cretion disc dShakura & Sunvaevll 1973b . and it follows closely the 
method that we introduced in Paper I, though we will now consider 
a general opacity law and equation of state. For this reason we only 
summarise the steps leading up to Eq. (21) of Paper I, after which 
we deviate from the treatment in Paper I. 

The surface density of the accretion disc is defined as 



-£ 



pdz = 2pH, 



(1) 



where p is the density, and H the half thickness of the disc. We can 
then write the accretion rate as 



M = -2nRLv R , 



(2) 



which is independent of the radial coordinate R, and where v R is the 
radial velocity. 

A thin accretion disc is supported by the centrifugal force in 
the radial direction, so we get that 

, GM 

4 ~ — = 0, (3) 

where is the azimuthal velocity, G the gravitational constant and 
M is the mass of the accretor. This shows that the disc rotates in a 
Keplerian fashion. 

The disc is in hydrostatic equilibrium in the vertical direction, 
which gives us that the pressure at the midplane of the disc is 



1 GM GMH 2 
P= -HI.— = p — 

2 R 3 F R 3 

but the hydrostatic equilibrium can also be expressed as 
H _ c s 

where 



(4) 



(5) 



c s 



1/2 



(6) 



is the isothermal speed of sound. This shows that the Keplerian 
velocity is highly supersonic in a thin accretion disc. 

The conservation of angular momentum leads to the following 
height-integrated equation (Paper I) 



d/ 
VR dR~ 



B-Bj, 



Ho 
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R+ 

RdR 



, d 
R 3 vl, — — 
dR\R 2 



(7) 



where the specific angular momentum / = Rv^, oc R 1/2 and v is 
the viscosity. The magnetic term describes the exchange of angular 
momentum between the disc and the star via the magnetosphere. 
We assume that the vertical magnetic field is due to the dipolar 
field of the neutron star, so that its value in the equatorial plane is 



B 7 = --£r 



R 3 



(8) 



where p is the magnetic dipole moment of the star. There are two 
sources for S^. The shear between the Keplerian disc and the mag- 
netosphere produces a field whose value in the upper half of the 
disc is 



£<4, S hear = ~jB z 1 - — 



(9) 



where fl k = v^/R, Q s is the angular velocity of t he star, and y is 
a dim ensionless parameter of the order of a few dGhosh & Lambl 
1 1979b . 

The magnetohydrodynamical turbulence in the accretion disc 
is also generating a magnetic field, B Aim , th rough a dynamo ac- 
tion (see for instance lBalbus & Hawlevll 1998b . In order to estimate 
the strength of this magnetic field we start by making the assump- 
tions that the turbulence is not affected by B^sheao which is justi- 
fied b ecause in general |5^, tU rbl > IS^shear] dTessema & TorkelssorJ 
2010), and that the accretion is driven in its entirety by the 
small-scale turbulent magnetic (Maxwell) stress. The latter as- 
sumption ignores the fact that the Reynolds stress is responsi- 
ble for roughly 20% of t he radial transport of angular momentum 
(Brandenbu rg et aljl995h , but this results in only a minor quantita- 
tive change. It is then the ^-component of the turbulent Maxwell 
stress tensor, (Bi,turb-B#,turb)/A'0! that transports the angular momen- 
tum, and we couple this t o the pressure in the disc through the 
Shakura & Sunvaevl Jl973b a ss -prescription 



Ho 



(10) 



Note that a negative Maxwell stress means that the angular momen- 
tum is transported o utwards and that numerical simulations su: 
gest that or ; s ~ 10 ~ 2 dHawlev et al.lll99llBrandenburg et aD 1993 
iTorkelssort (l998) suggested thaQ 
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where y dvn ~ 1 bas ed on the numerical simulations by 
Brandenb u7g et all dl995h . We now assume that due to an in- 
verse cascade this turbulent filed contains a large-scale field, B^ yn 
dFrisch. Pouquet. Leorat & Mazurelll975h . and we express the ratio 
of the large-scale field to the rms-value of the turbulent field as 



1 This is a more pre cise definition of yi yn than = yi yn B, that was used 
bv lTorkelssoi] d 19981) . 



Thin accretion discs around millisecond X-ray pulsars 3 



.1/2- 



(12) 



0,turb/ 



Thus we can express the large-scale magnetic field, that couples to 
the magnetic field of the neutron star, as 

/ \l/2 

B#,ayn = e\a ss noydynP{r)) , (13) 

where -1 < e < 1, and a negative value describes a magnetic 
field which is pointing in the negative (^-direction at the upper disc 
surface. We will in general be conservative and assu me that |e| < 
0.1, but the simulations by Brandenburg et "all i 19951) suggest that 
|e| can be significantly larger. 

We can now re-express Eq. Q as 
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(14) 



which can be formulated as an ordinary differential equation in vS 
(Paper I). 

2.2 Heating and radiative transport 

For a slow inflow of matter through an optically thick disc the local 
viscous dissipation v ■ f„ is balanced by the radiative losses V • F ra d. 
This leads to the height-integrated equation 

4crr L 4 9 GM 

= -vS— , (15) 

3t S R 3 

where T c is the temperature at the midplane of the disc, <x the 
Stefan-Boltzmann constant, and the optical depth of the disc is 
given by: 



Jo 



(16) 



where k is the sum of the electron scattering opacity 

k cs = 0.04 m 2 kg-', (17) 

and the free-free opacity, which is given by Kramer's law 

% = / f0 pr c - 7/2 m 2 kg- 1 (18) 
with 

k = 5 x 10 20 m 5 kg- 2 K 7/2 . (19) 
Combining Eqs. (14) and (15), we get 

d 27 GM 
Tt = —pH(K R + K es )vX— (20) 
32<x 



2.3 The vertical structure of an accretion disc 

The total pressure is the sum of gas and radiation pressure 

P(P,T C ) = + — -1 (21) 

m py u 3c 

where k B is Boltzmann's constant, fi = 0.62 the mean molecular 
weight for a fully ionised gas, m p the mass of a proton, and c the 
speed of light, but the pressure can also be expressed using Eq. l[4j 
for hydrostatic equilibrium. We can thus express the scale height in 
terms of the total pressure 

r4 P 3\i/2 

(22) 



[mpfiGM + 3cpGM) 



The viscous stress tensor (Paper 1, Eq. 30) gives us 
a ss P, 



3 gv /GM \ 
4~H\1F) 



(23) 



but pressure and density are related according to Eq. so we can 
solve for the density of the gas 



3 vS (GMy [/2 

Combining Eqs. 10, d!4t and &241 and using 
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where 



GMP 



4n 2 



1/3 



1.5 x 10 6 P 2/3 M, 1/3 1 

spin 1 



(24) 
vL we get 

(25) 
(26) 



is the co-rotation radius, at which the Keplerian angular velocity is 
the same as the stellar angular velocity. Here P sp i„ = |p is the spin 
period of the star and Mi = . As R — > oo 

M 
3n 

which we use as the outer boundary condition. 

We now introduce the dimensionless variable A through 

y = AM, (27) 

and a dimensionless radial coordinate through 

R = rR A , (28) 

where R A is the Alfven radius, which is giv en by setting the mag- 
netic pressure equal to the ram pressure (e.g. lShapiro & Teukolskvl 
Il983h 



( 2nY 

[gmm 2 ^ 



1.4 x WM\l"M^ 



(29) 



where M 14 represents the mass transfer rate in units of 10 14 kg s _1 , 
and yu 16 is the stellar magnetic moment in units of 10 16 Tm 3 . The 
Alfven radius should not be taken as the location of the inner edge 
of our discs, since this is a quantity that comes out of our solutions, 
but merely as a convenient length unit. Furthermore we introduce 
the fastness parameter 



: 0.36M- 5 ' 7 M- 3 ^^ 



4.8 ms 



(30) 



We can now write Eq. l !25t ; 



dA 
~dr 
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6nr 
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^-1/4^-5/4^-9/4 



R -V2 r -W 



[l- Ws r 3 ' 2 ], (31) 



3 REGIONAL DISC STRUCTURE EQUATIONS 

In order to solve Eq. (fJTJ numerically we need to determine H for 
the different regions in the disc. With gas pressure and Kramer's 
opacity we get 
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and for the gas pressure and electron scattering opacity 

/ 81*^ V /10 / k 3 \ 2/5 /(r# A ) 3 \ 7/2 \ ■ W 

H(r) = — — — — MA . (33) 

\mo-a s J \m p /X) \ GM J V J 

Finally with radiation pressure and electron scattering the scale 
height is 

g = -Jlyj, = — ^(MA) (34) 

oC oC 

which is essentially independent of R. This comes from the fact 
that for radiation pressure T* ~ v2. 2 MR~ 3 . 

We can now write down the self-similar structure of the differ- 
ent regions of the accretion disc in terms of the function A(r). The 
toroidal magnetic field is given by the same expression 

iW = -3.6 x lOVM^M^V^ 7 (l - ^n-) r ~ 3 T (35) 

everywhere in the disc, but for all other quantities we have to handle 
the different regions of the disc separately. 

In the outer disc where P g » P ra d, and Kg s> k cs we have that: 

I = 1.6 x 10V /4 ^ /5 Mf 4 M^ /3 Vi 6 3/7 A(r) 7/10 r- 3/4 kg nT 2 (36) 

r c = 8.1 x io 7 /; 1/ V /5 ^ /14 ^ 1 '4 8/3 V7 6 3/7 a(^) 3/10 ^ 3/4 k (37) 



v R = 1.6 x 10 3 /i- 4/5 a 4 / 5 M7 1/7 M; 4 8/3 Vr 6 8/35 A(r)- 3/5 r- 2/5 m s" 1 (50) 



H 



- = 8.1 x 10- 3 /i- 3/8 < s 1/,0 M7 11/28 M 4 / 3 V!r A ('') 3/20 - 1/8 



(38) 



Pc = 7.1 x 10 3 // 9/8 < s 7/10 M 2S/28 M 3 4 8/3 V7 6 15/14 

A(r )ll/20 r -15/8 kgm 3 



(39) 

M (40) 

y = 6.4 x 10 7 /i- 3/4 a 4/5 M- 5/14 M 3 / 3 V 3 ^ 7 A(r) 3/1() r 3/4 m 2 s" 1 (41) 



r = 5.9 x 10 2 /jff- 4/5 M.'f A(r) 1/5 



v R = 1 .4 x 10 2 /T 3, V 4 f MT 3/14 Mjf 5 //- 1/7 

A(r) -7/10 r -l/4 ms -l 



zw = 7.7 x lov/^r ^r^r 

A(r) 17/40 r -21/16 T 



(42) 



(43) 



The transition from the outer to the middle region occurs where 
= Kg. Assuming that A = 1/3/r at this point we have that the 
transition radius is 

,1-7 — 1/3 J^10/21 »V20/21 "4/7 

which is independent of a ss . 

In the middle disc, where P g » P mi , and k cs » we have 

that: 

I = 7.1 X ^/^a^^Mf ^ 7/3 V7 6 12 ' 35 A(r) 3/5 ^ 3/5 kgm- 2 (44) 



T c = 1.8 x 10 /< ' a~s ' M 3 ' M 23 ' ' A(r) ' r ' K (45) 
| = 1.2 x 10- 2 /2- 2/5 «-' s 1/1 X 5,14 *i4 /7 V|f A(r)'/V 1 / 20 (46) 



Pc = 2.o x io 3 /} 6/5 < 7/10 m; 1/I4 m« 4 /7 V7 6 33/35 

A(r) 2/5 r -33/20 kgm -3 

t= 1.4 x 10 4 /i 13/1 V /5 < 7 * 2 4 7/3 V 3 fA('-) 3/5 ^ 3/5 
v = 1.5 x 10 8 /2- 4/5 a 4 / 5 M- 2/7 M 8 / 3 V|6 /35 A(r) 2/5 r 3/5 m 2 s" 1 (49) 



(47) 
(48) 



fi ,, dyn = 6.3 x io 4 ^'/ 5 ^^ 20 ^; 7 ' 28 ^;: 7 ' 14 ^^ 1770 



A(r) 2/5 r -5 1 /40 T 



(51) 



The disc structure changes drastically in the inner disc where 
Prad s> Pg and K es » Kg, where we have that 

I = 9.6 x 10V; s 1 M 1 " 5/7 M M 10/ Vj 77 A(r)-'r 3/2 kgm- 2 (52) 



T c = 1.9 x 10 7 ^ /4 Mf 28 M 3 ' 2 V 6 3/14 r- 3/8 K 



| = l.l«^A(r)r-' 



(53) 
(54) 



p c = 3.2 x 10- 3 a ! : s 1 M7 5/7 M7 4 17/ V* 77 A(f-)- 2 r 3/2 kg nT 3 (55) 
T„ = 1.9^ 1 M- 5/7 M- 4 10/ V^ 7 A(r)- 1 r 3 / 2 (56) 
y = 1.0 x 10 I2 a 8S Mf 7 M; 4 7/ V7 6 6/7 A(r) 2 r- 3/2 m 2 s" 1 (57) 



(58) 



v R = 1.2 x 10 7 a ss Mf 7 M 1 ' 4 9/ V7 6 10/7 A(r)r- 5/2 ms- 1 

S,, dyn = 6.6 x lOWjyrMl' 14 M 3 ^V: 3/7 r- 3/4 T (59) 

It is of particular significance here that S oc A -1 . 

The transition radius between the middle and the inner region 
is estimated by equating the gas and radiation pressures in the mid- 
dle region and approximating A with 1/3tt 

r IM = 12^ ^M^Mf^^J 1 

We then see that the accretion disc right outside of the Alfven radius 
is dominated by radiation pressure only if 



(60) 



This condition is not fulfilled for a conventi onal X-ray pulsar wit h 
a magnetic dipole moment of ~ 10 20 T m 3 dWhite & Stelialll988h . 
though it can be fulfilled for a millisecond X-ray pulsar. 



4 GLOBAL SOLUTIONS 

For our fiducial model we take a neutron star of M = 1 .4M Q with a 
magnetic moment of 10 16 Tm 3 and a spin period of 4.8 ms. We set 
the dimensionless parameters y and y dyn to 1 and 10, respectively, 
while a S!i = 0.01. y, ydyn. and a ss appear only in combinations with 
other parameters, and thus their values do not carry any special 
significance, but the parameters e and M appear on their own in 
the equation, and therefore influence the solution in unique ways. 
In particular M influences the fastness parameter and the transition 
radii between the regions of the accretion disc. In this paper we 
consider three different accretion rates, M l4 = 0.012, 0.12 and 1.5. 
We list the corresponding fastness parameters and transition radii 
in Tab. □ 



4.1 Case I:R A > ^om 

The accretion disc consists of only an outer region when r 0M < 1, 
which corresponds to that M < 1.4 x 10 I2 kgs~' for our fiducial 
model. Assuming M = 0.012M 14 we get the solutions from the top 
to the bottom of Fig.[[]for e = 0.1, 0.05, 0, -0.05 and -0.1. We 
see here that A — > at a finite r for e = -0.1, and -0.05, which 
we called case D in Paper I, since p — > at this radius (Fig.[2j. In 
the other cases A grows without a bound for decreasing r, which 
we call case V in Paper I, since then the viscous flux of angular 
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Table 1. The Alfven radius, fastness parameter a> s and the transition radii in 
units of the Alfven radius for different accretion rates M14 and the fiducial 
neutron star. 



M14 R A [m] oj s r M mi 



0.012 


4.7 X 10 4 


0.98 






0.12 


2.4 x 10 4 


0.37 


8.5 




0.8 


1.4 x 10 4 


0.31 


52 


6.0 


1.5 


1.2 x 10 4 


0.12 


95 


12 



0.20 



0.15 - 



< 0.10 - 



0.05 - 




0.00 1 

10 20 30 40 50 



Figure 1. A(r) for our fiducial neutron star with an accretion rate of 1.2 X 
10 12 kgs-', and e = 0.1, 0.05, 0, -0.05, -0.1 from the top to the bottom. 
In this case the disc consists of only an outer region that is dominated by 
gas pressure and Kramer's opacity. 

momentum changes sign at a finite radius where j-(-y'rA) = 
(Fig.O, which we take as the inner edge of the disc. Inside of this 
radius the accretion flow is driven by the magnetic coupling to the 
neutron star. 

4.2 Case II: R 1M < Ra < Rom 

The disc has both an outer and a middle region when 0.014 < 
M\4 < 0.14. We then use the analytical solution with A = l/3n 
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Figure 2. Z as a function of r for our fiducial neutron star with an accretion 
rate of 1.2 x 10 12 kgs~', and with e = 0.1, 0.05, 0, and -0.1 from the top 
to the bottom. 



0.10 rr 
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5 10 15 20 

r 

Figure 3. ^ VPa) as a function of r for the fiducial neutron star with an 
accretion rate of 1.2 X 10 12 kg s , and with e = 0.1, 0.05, and from the 
bottom to the top. 
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Figure 4. A(r) for the middle region of an accretion disc around our 
fiducial neutron star with an accretion rate of 1.2 X 10 13 kgs~', and with 
e = 0.1, 0.05, 0, -0.05, -0.1 from the top to the bottom. 

for the outer region and solve Eq. t3U numerically for the mid- 
dle region. For our calculations we take M = 1.2 x 10 13 kgs~', 
which yields tom ~ 8.5. The numerical solutions for the mid- 
dle region are shown in Fig. [4] for our fiducial neutron star with 
6 = 0.1, 0.05, 0, -0.05, and -0.1. We find a case D inner bound- 
ary for e = -OA, while the other solutions have case V inner bound- 
aries. The transition from case D to case V at the e, that gives the 
smallest possible radius for the inner edge of the disc, which varies 
weakly with the accretion rate M. 

4.3 Case III: R A < R m 

The innermost region of the disc is dominated by radiation pressure 
and electron scattering when M14 > 0. 14. We use two different ac- 
cretion rates that fulfill this condition, M = 8 x 10' 3 kgs~' and 
M = 1.5 x 10 14 kg s~' together with our fiducial neutron star in our 
calculations. We assume that A = 1/3^ in the outer and middle 
regions of the disc, and solve Eq. d3U for the inner region starting 
from Tim s; 12. Our solutions for the inner disc region at the higher 
accretion rate and with e = 0.1, 0.05, 0, -0.05, -0.1, and -0.12 
are shown in Fig. [5] All these solutions have case V inner bound- 
aries except the e = -0.12 model. However because £ oc A~' in 
this case we see that the density goes to infinity at the inner edge of 
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Figure 8. ( VrAj as a function of r for the fiducial neutron star with an 
Figure 5. A(r) for our fiducial neutron star with an accretion rate of 1.5 X accretion rate of 1.5 X 10 14 kg s~\ and with e = 0.1, 0.05, 0, -0.05, and 

10 14 kgs~', and with e = 0.1, 0.05, 0, -0.05, -0.1, -0.12 from the top to -0.1 from the bottom to the top. 

the bottom. 



2.0x10" - 




r 



Figure 6. £ as a function of r for the fiducial neutron star with an accretion 
rate of 1.5 X 10 14 kgs"', and with e = 0.1, 0, -0.1, and -0.12 from the 
bottom to the top. 
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r 

Figure 7. Vr as a function of r for the fiducial neutron star with an accretion 
rate of 1.5 X 10 14 kg s" 1 and e = 0.1, 0, -0.1 and -0. 1 2 from the top to the 
bottom 



the e = -0.12 disc (Fig. [6), at which point the flow stagnates (Fig. 
[7}. In the case V solutions, on the other hand, we see that X — > at 
a point which is inside the inner edge of the disc (Fig. [8}. 



5 DISCUSSION 

5.1 The inner edge of the accretion disc 

In order for the neutron star to appear as an X-ray pulsar the accre- 
tion disc must be truncated by the stellar magnetic field (Tab.[2} at a 
radius R outside of the neutron star. The inner edge of the disc can 
be either of type D (A = 0) or type V ( ^-( -\JrA) = 0) (see Paper I 
for an extensive discuss ion of these bou ndaries). In models without 
a dynamo (for instance [Ghosh & Lamb|[l979l) Ro is comparable to 
or smaller than the Alfven radius, ^a, but we find that the dynamo 
can make it significantly larger than R A . This opens up the possibil- 
ity, as suggested by the referee, that a system can temporarily turn 
into an X-ray pulsar during the periods when the magnetic field 
generated by the dynamo is strong enough to make Rq larger than 
the radius of the neutron star. S AX J1748. 9-2021 is suc h an inter- 
mittent millisecond X-ray pulsar jAltamirano et al . 2008), which is 
on for hundreds of seconds at a time. 

In case R is located in the outer or middle region of our disc 
solution we find that a type D boundary coincides with the point at 
which p — > and v R — > oo though the accretion rate is finite. How- 
ever a type D boundary in the inner part of the accretion disc be- 
haves in the opposite way; p — > oo and Vr — > 0. Neither one of these 
behaviours is entirely physical and one must keep in mind that we 
have ignored the importance of the radial p ressure gradient and the 
inertial term in Eq. Q (s ee for instance [Abramowicz et al.ll 19881 : 
Poph am & Naravarj Il992l) . and we have not attempted to model 
how the plasma is entering the stellar magnetic field lines to form a 
funnel flow onto the magnetic poles of the neutron star. One further 
problem with the inner region of the di sc solution is that it is know n 
to be viscously and thermally unstable (Lig htman & Eardlevl Il974l) . 
This instability pe rsists if the accre tion disc is threaded by a stel- 
lar magnetic field dCampbell|[l998h . but so far it has not been in- 
vestigated what happens if the accretion disc is sustaining its own 
magnetic field through a dynamo. 
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5.2 Accretion torque 

To understand the nature of the exchange of angular momentum 
between the accretion disc and the neutron star we multiply Eq. 
d!4t by 2nR and integrate it from R , the inner radius of the disc, to 
a large radius Ri. 



-M^yjGMRt - ^GMR ) = 

JR A<0 Jr A<o 

-3n(vS) Rl ^GMR 1 + 3n(vT) Ro ^GMR Q . 



(61) 



The left-hand-side is the difference between the angular momentum 
that is adverted through the inner edge of the accretion disc and that 
which is fed into the disc at its outer edge and the right-hand-side 
describes the contribution of magnetic and viscous torques to the 
angular momentum balance. Notice that the term 3n(vS) Rl s/GMRi 
describes the viscous tension at the outer edge of the disc, which 
will not be considered further in this paper. The material torque of 
the inner edge of the disc on the neutron star is 



(62) 



N adr = M^GMR = 1.4 x \0 26 ^ M 3 / 1 M% 7 r^' 2 . 

It increases as the accretion rate increases (TabO. Then we have 
the magnetic torques on the neutron star, which we have divided 
into one part due to the shear 



B-Btp shear n 

An R~dR ; 



r 

3.6 x 10 2 VM /7 M^ 
r[r- 4 (l-c s r 3 ' 2 )]dr 

and a second part due to the dynamo 



N, 



dyn 



JRn 



An 

Ro A*0 



dyn „2 



R z dR. 



(63) 



(64) 



We split up the latter integral into one integral for each of the three 
different regions of the disc 

pouter = 7.7 X W^eyfa^^^M^M^ 



t-'OUt 



A 17/4 V 37/16 rfr, 



N, 



dyn.mid ' 



,27^,1/2^1/20-1/5,, 19/70 ^17/28^107/140 



16 K 



14 



A 2/5 r- 91/40 dr, 



(65) 



(66) 



and 



N, 



dyn,inner 



6.6 x 10- 



^^MrMT f r-^d, (67) 

*/ inner 



Finally we have the boundary terms describing the viscous stresses 
on the inner and the outer edges of the accretion disc. The viscous 
stress at the outer edge of the disc is responsible for removing the 
angular momentum that is transported outwards through the accre- 
tion disc, and does not affect the angular momentum evolution of 
the neutron star. The viscous stress at the inner edge of the disc on 
the other hand can result in a torque that transports angular momen- 
tum away from the neutron star 



N vm = -3n(vI.) Ro y[GMR = 
-1.29xl0 2 y 6 7 Mf 7 M^A(,- o) ^ 2 



(68) 



This term is usually neglected in accretion disc theory since the 
standard accretion disc solution has a case D inner boundary. But it 
can transport angular momentum from the neutron star to the disc 
when we have a case V inner boundary. 

We report on all the components and the total sum of the 
torques on our fiducial neutron star in Tab[2] The main factor de- 
termining the total torque is the accretion rate. In Paper I we found 
that the dominant contributor to the torque is the internal disc dy- 
namo. This torque is still strong, but now A? adv and N vis are compa- 
rable in strength, though these make opposing contributions to the 
net torque, while the contribution from is still small. There is 
therefore no simple dependence of N tota j on e, though the dynamo 
is indirectly affecting /V a dv and /V v i s by causing the inner d isc edge 
to move outwards compared to the iGhosh & Lambl d 19791) model. 
/Vyis is also highly sensitive to the nature of the inner disc boundary, 
and in particular it vanishes for a D inner boundary, which occurs at 
a sufficiently negative e, in which case the net torque is determined 
by the balance between and N^ y . 

5.3 Observed properties of millisecond X-ray Pulsars 

The largest spin up rate recorded so far for a millisecond X- 
ray pulsar is ~ 10~ 12 Hz s" 1 during the D ecember 2004 outburst 
of IGR J0029 1+5934 jBurderi et all l2007h . but one should keep 
in mind that there are large uncertainties in the spin variations 
that hav e been reported for the millisecond X-ray pulsars. For 
instance iBurderi et all d2006t) reported vs between -7.6 x 10~ 14 
and 4 4 x 10 -13 Hzs -1 for SAX J1808.4-3658, but lHartman et all 
(2008) noted that the measurements of this source are plagued 
by large variations in the pulse shape, and put an upper limit of 
2.5 Xl0 H Hz s on the spin variations at any instant. Some of the 
timing noise in the millisecond X-ray pulsars can be explained b y 
the motion of the hot spot on the neutron star jPatruno et al.ll2009l) . 
and may thus not correspond to real changes in the spin period of 
the neutron star. 

The spin variation corresponds to a torque 



N = 2nvl = 6.3 x 10 23 v 13 /3 8 Nm, 



(69) 



where I$g is the moment of inertia of the neutron star measured in 
10 38 k gm 2 , and Vp is the spin derivative in units of 10~ 13 Hzs~'. 
In the IGhosh & L amb ( 1979) model, the spin variation of IGR 
J0029 1+5934 would imply then an accretion rate of at least ~ 
10 14 kg s~', which is uncomfortably large, but our Tab.|2]shows that 
it is possible to enhance the torque severalfold by coupling the stel- 
lar magnetic field to the magnetic field generated by an internal 
dynamo in the disc, and through the fact that the inner edge of the 
accretion disc is then located outside the Alfven radius, which in- 
creases the amount of angular momentum that is adverted onto the 
neutron star. Thus it is possible that a lower accretion rate can pro- 
duce the necessary torque, though we have not attempted to fit our 
model to this system. 

5.4 The validity of the model 

A significant part of the luminosity is not generated in the accretion 
disc outside of R , but rather it is released as the accreting matter 
falls onto the neutron star. We estimate that this luminosity is 



: GMM 



1 

Rut 



1 



1 GMM 



2 R 



GMM 



1 



2R ) 



(70) 



where we have assumed that both the potential energy difference 
between the inner edge of the accretion disc and the surface of the 
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Table 2. The inner edge of the accretion disc and its torque on the fiducial neutron star for different accretion rate. 



M14 


e 


Rolrn] 


A'dyn.ouler [Nm] 


Wd yn ,Middle [Nm] 


^dynanner [Nm] 


N sh [Nm] 


Wadv [Nm] 


N V1S [Nm] 


Metal [Nm] 


0.012 


0.1 


1.8 X 10 5 


3.4 x 10 24 


- 


- 


-7.8 x 10 23 


7.0 x 10 24 


-7.4 x 10 24 


2.2 x 10 24 




0.05 


9.9 x 10 4 


3.7 x 10 24 


- 


- 


-1.7 x 10 24 


5.2 X 10 24 


-6.7 x 10 24 


5.0 x 10 23 







4.7 x 10 4 





- 


- 


-3.1 x 10 24 


3.6 x 10 24 


-3.4 x 10 22 


4.7 x 10 23 




-0.05 


8.8 x 10 4 


-2.9 x 10 24 


- 


- 


-1.9 X 10 24 


4.9 x 10 24 





1.0 x 10 23 




-0.1 


1.2 x 10 5 


-4.1 x 10 24 


- 


- 


-1.4 x 10 24 


5.6 x 10 24 





1.0 x 10 23 


0.12 


0. 1 


9.4 X 10 


8.4 X 10 


1.4 X 10 




-2.1 x 10 


5.1 x 10 


-l.O X 10 


4.7 X 10 




0.05 


6.0 x 10 4 


A 1 v f 1 Ail 

4.1 X 10 


1.5 x 10 25 




-3.2 X 10 


4.1 X 10 


-7.7 X 10 25 


1.7 X 10 25 







3.0 x 10 4 










-7.4 X 10 


2.9 X 10 


-2.4 X 10 


4.3 X 10 




—0.05 


3.2 x 10 4 


—4.0 X 10 


-2.6 x 10 25 




—2.0 X 10 


T n , j i ^25 

2.9 X 10 


—6.7 x 10 22 


o f\ , j 1 a25 
—3.9 X 10 




-0. 1 


4.8 X 10 


-7.9 X 10 


-3.0 x 10 25 




—3.7 X 10 


3.6 X 10 25 





-7.7 X 10 25 


0.8 


0.1 


3.2 X 10 4 


4.9 X 10 24 


3.1 x 10 23 


8.2 X 10 25 


-1.0 x 10 25 


2.0 X 10 26 


-2.7 x 10 26 


7.2 X 10 24 




0.05 


1.6 X 10 4 


2.5 x 10 24 


3.2 x 10 25 


1.1 X 10 26 


2.7 x 10 25 


1.4 x 10 26 


-2.5 x 10 26 


6.2 x 10 25 







1.7 X 10 4 











1.4 x 10 25 


1.5 x 10 26 


-1.6 x 10 26 


4.0 x 10 24 




-0.05 


2.0 X 10 4 


-2.4 x 10 24 


8.3 x 10 22 


-7.7 x 10 25 


2.5 x 10 23 


1.6 x 10 26 


-6.2 x 10 25 


1.9 x 10 25 




-0.1 


1.7 X 10 4 


-4.8 x 10 24 


-1.3 x 10 25 


-1.8 x 10 26 


1.4 x 10 25 


1.4 x 10 26 





-4.4 x 10 25 




-0.12 


1.8 x 10 4 


-5.8 x 10 24 


-8.4 x 10 24 


-2.0 x 10 26 


-3.2 x 10 24 


1.5 x 10 26 





-6.7 x 10 25 


1.5 


0.1 


2.6 x 10 4 


2.1 x 10 24 


2.5 x 10 25 


1.1 x 10 26 


-5.6 x 10 24 


3.3 x 10 26 


-4.5 x 10 26 


1.2 x 10 25 




0.05 


2.1 x 10 4 


1.0 x 10 24 


1.3 x 10 25 


7.1 x 10 25 


4.6 x 10 24 


2.9 x 10 26 


-4.2 x 10 26 


-4.0 x 10 25 







1.8 x 10 4 











1.7 x 10 25 


2.7 x 10 26 


-3.2 x 10 26 


-3.3 x 10 25 




-0.05 


1.6 x 10 4 


-1.1 x 10 24 


-1.3 x 10 25 


-9.1 x 10 25 


3.7 x 10 25 


2.6 x 10 26 


-3.2 x 10 26 


-1.3 x 10 26 




-0.1 


1.5 x 10 4 


-2.1 x 10 24 


-2.5 x 10 25 


-2.0 x 10 26 


5.5 x 10 25 


2.5 x 10 26 


-1.1 x 10 25 


6.7 x 10 25 




-0.12 


1.3 x 10 4 


-2.5 x 10 24 


-3.0 x 10 25 


-2.4 x 10 26 


5.8 x 10 25 


2.5 x 10 26 





3.6 x 10 25 



neutron star as well as the kinetic energy associated with the Kep- 
lerian motion at the inner edge of the accretion disc is released at 
the surface of the neutron star. Consequently GMM/2R NS < L NS < 
GMM/R m . 

The neutron star can then contribute to the local heating in the 
accretion disc through its irradiation of the disc. Assuming that the 
neutron star is a point source and the accretion disc is flaring, i.e. 
din H/dln R > 0, the irradiation flux on the accretion disc is at 
most dFrank. King & Rainell2002h 

L NS H /din// \ 

-1(1-/3) = 



(71) 



4^ 2 R \ d\riR 
GMM HldlnH 



4nR m R 2 R\d\nR 



-1 (1-/3), 



where {3 is the disc albedo (for a more complete treatment taking 
into account shadowing by the disc see Ed ams &ShJl98fj) . The 
local flux generated by the viscous heating in the accretion disc is 
on the other hand 



_ 9 GMVL 
Faccr " 8~^' 
We thus find that 



2 H /dln# 
3~R^ \ dlnfl 



1 (1-/3), 



(72) 



(73) 



and since dln/f/dlni? - 1 ~ 0.1 we find that irradiation dominates 
for H > 15Z?ns if j8 = 0. Consequently irradiation dominates over 
viscous heating for R > 10 4 km. Thus the outer part of the accretion 
disc is hotter than is pred icted by our model, which can be observed 
as an increased UV flux JVrtilek et alF l990). Thus we have under- 
est imated the pressure i n the outer disc, and since we have used 
the lShakura & Sunvaevl ll973l) prescription ct^P to estimate S^dyn, 
we have also underestimated B^dyn and the torque from the outer 
disc on the neutron star. However, this torque is small in compari- 
son to the total torque on the neutron star since S-B^ dyn is a rapidly 



decreasing function of R, and thus the effect of the irradiation of 
the outer disc does not significantly change our results. Irradiation 
is unlikely to significantly affect the inner disc, but, as pointed out 
above, other effects, such as radial advection and the transition from 
a disc flow to a funnel flow, that are not considered in this paper be- 
come important in the inner disc. We will investigate these effects 
in the future. 



6 CONCLUSIONS 

The accretion rates that have been observed in millisecond X-ray 
pulsars cover several orders of magnitude, and the highest accretion 
rates are sufficient that the innermost part of the accretion disc is 
dominated by radiation pressure and electron scattering. For these 
reason s we have followed the approach by IShakura & Sunvaevl 
dl973l) and divided our disc model in three regions, an outer re- 
gion dominated by gas pressures and Kramer's opacity, a middle 
region dominated by gas pressure and electron scattering, and an 
inner region dominated by radiation pressure and electron scatter- 
ing. 

The large spin changes that have been observed in some of 
the accreting mi llisecond X-ray puls ars are difficult to explain in 
the model by lGhosh&Lambl i ll979h , but the inclusion of a disc 
dynamo produces a significant enhancement of the torque between 
the accretion disc and the neutron star. This increase is not only due 
to the coupling between the magnetic fields of the neutron star and 
the accretion disc as such, but also to the fact that the accretion disc 
is truncated at a larger radius thus increasing the lever arm between 
the edge of the disc and the neutron star. 
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